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A method for estimation of the heat-shielding properties of gas-fdled mineral wool insulators is considered, 
which takes into account the mechanism of interrelated heat and mass transfer processes and their intensities, 
depending on the service conditions, design specifics, and technological parameters of the insulating material. 


Heat exchange in fibrous heat-insulating materials which 
constitute highly disperse gas-filled systems, is a complex 
process which includes heat transfer by conduction, radia¬ 
tion, and convection and is implemented by means of: 

- conductive heat transfer in the fiber which forms the 
matrix of the porous structure of a heat-insulating material; 

- radiative heat exchange inside the absorbing, reflect¬ 
ing, and dissipating medium of the inter-fiber space; 

- convective flows of the gas filling the inter-fiber space, 
caused by the temperate difference between the boundary 
surfaces of the heat-insulating structure. 

The heat-shielding properties of fibrous materials are 
usually determined by employing relationships based on the 
Fourier law, which, strictly speaking, are only true for homo¬ 
geneous solid bodies in which heat transfer in the presence of 
a temperature gradient proceeds only in the form of heat con¬ 
duction. 

It is obvious that heat transfer in such an approach is a 
generalized characteristic of the thermal conductivity of the 
material, i.e., the effective thermal conductivity X f : 

Kf=\ + K + h’ C 1 ) 

where X , X r , and X c are conductive, radiative, and convec¬ 
tive heat transfer, respectively. 

The present study, on the basis of the effective heat trans¬ 
fer concept, considers a method for the estimation of heat¬ 
shielding properties of gas-filled fiber insulation, which is 
based on contemporary concepts and takes into account the 
mechanism of interrelated heat- and mass-exchange pro- 
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cesses and their intensities, depending on the operating con¬ 
ditions, structural specifics, and technological parameters of 
the insulating material. 

Conductive heat transfer. An analysis of data obtained 
on the basis of the variation principle of the continuous me¬ 
dium mechanics and other theoretical models for conductive 
heat transfer, as well as semi-empirical models including ex¬ 
perimentally determined coefficients, demonstrated that in 
using these coefficients, one can obtain values which are the 
closest to the experimentally found values of thermal con¬ 
ductivity [1-5], 

Let us consider such a model described in [5], which 
gives satisfactory results in the calculation of conductive 
heat transfer in heat-insulating materials extensively used for 
heat-insulating structures in industrial insulation. 

The calculation formulae of this model used for determi¬ 
nation of the conductive heat transfer for insulation in the dry 
state have the following form: 

X,, =---; (2) 

7 (1 - a)/X l + a/X 2 

1 (1 -m)X t + mX g 

= H ^ n W 

(1 -m)X f + m /X g 

where a is the structural coefficient (found experimentally); 
m is the porosity of the insulating material, fraction of unity; 


314 


0361-7610/00/0910-0314$25.00 © 2001 Plenum Publishing Corporation 




A Method for Estimating Heat-Shielding Properties of Mineral Wool Insulation 


315 


A, f is the heat conductivity of the fiber; is the heat conduc¬ 
tivity of the gas filling the porous structure of the material. 

The proposed model also makes it possible to determine 
the conductive heat transfer of a moist material: 


q (1 -a^ + a/X”’ 


1 (l-m)A, f + WX f +{m-W)X 2 

__ [ _ 

2 ~~ (l-/w)/^- av + W /X w + (m -W)/ Xg ’ 

where X is the heat conductivity of water; W is the relative 
volume moisture of the insulation related to the relative 
moisture content of the material U (kg of moisture per kg of 
dry mass) by the relationship W= f/y/p w (p w is the water 
density; y is the volume weight). 

Based on formulae (2) - (4), it is possible to calculate the 
conductive heat transfer at any temperature, taking the values 
A, f and X at the average temperature of the insulation t, w . 

The conductive heat conductivity of a moist material can 
be determined subject to the condition: 
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Fig. 1 . Dependence on temperature of values of heat transfer by ra¬ 
diation in fiberglass heat-insulating materials estimated according to 
methods [12] (1 ), [7] (2), [8] (3, 4 ), [9] (5), and [10, 11] (6). 


5°C < t m < 25°C. 

Heat transfer by radiation. An analysis of domestic 
and foreign publications [6 - 12] demonstrated that the rela¬ 
tionships describing the radiant heat conduction process in 
fibrous materials all have the same structure: 

K = f( d f. m ' a > b ’ c ) kT m ’ 

where d { is the fiber diameter; a, b, and c are the radiation 
properties of the fibrous medium; k is the Boltzman constant; 
T is the mean temperature of the insulator. 

For comparative analysis of the methods for determining 
the radiation heat transfer in fibrous materials proposed by 
various researchers, radiative thermal conductivity was cal¬ 
culated for fiberglass insulator with the volume weight 
50 kg/m 3 and the fiber diameter 6 x 10 ~ 6 m in accordance 
with the calculation relationships recommended in [7 - 12] 
(Fig. 1). A significant spread in the data is observed, nearly 
by the order of magnitude. The best convergence is seen in 
the results of methods [10 - 12]. 

This fact appears to indicate that the best approximation 
to the actual radiation heat transfer in fibrous insulating ma¬ 
terials is provided by the model suggested in [12]: 


where K at is the attenuation coefficient, which is a constant 


determined for various types of fibers by empirical formulae: 
mullite-silica fiber 

K at = 25mdjT m + 5 x 10 3 (1 — m ) 2 /T m ; 
basalt fiber 

K at = 70/W ( T m + 5x 10 3 (1 - m ) 2 / T m ; 

mineral wool fiber 

K at = 20md f T m + 7 x 10 3 (1 - m ) 2 /T m ; 
glass fiber 

K at = 11.2 x I06md f /T m + 800(1 - m ) 2 /T m . (6) 

The great practical significance of this model consists in 
the fact that it contains calculation relationships which take 
into account fiber properties. 

Heat transfer by convection is defined as heat trans¬ 
ferred by convective flows of the gas filling the porous struc¬ 
ture of the insulating material, which arise due to the differ¬ 
ence of temperatures on the boundary planes of the heat-in¬ 
sulating structure. 

A convenient method for studying these processes is the 
numerical solving of a system of nonlinear differential equa¬ 
tions describing the joint heat and mass transfer in porous 
media. 
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The system of equations describing the stationary con¬ 
vection process in a flat heat-insulating layer of height h and 
thickness 6 with impermeable lateral boundaries at which the 
temperatures 7j and T 2 are maintained (T 2 >T X ), using the 
Darcy linear filtration law 

V f = -grad p (7) 

d 


satisfied), the following criterion equation was found based 
on numerical modeling data: 

E k = (Ra® ) 0035 , (13) 

and for cylindrical heat-insulating structures [14]: 

E k = 1 + 5.9 x 10“ 3 [q(l — q)] 3/4 exp ((-4.29r|)(Ra®) 2 , (14) 


and approximation of the lift forces in the Bussinesque ap¬ 
proximation, has the form [13]: 


K 



P&PAT; 


( 8 ) 
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dp 

+ pg v pAr ; 

ay 


dU { 

dx 


dV f 

~dy 
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r/ dr" 
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d 2 T N 

U f — 

+ V f — 

| = A,® 

H - 7T“ 

t f dx 

dy y 


K dx 2 

dy ] 


(9) 

( 10 ) 


(ID 


where p is the pressure; p is the dynamic viscosity; K is the 
permeability coefficient; U { and V f are the filtration rate com¬ 
ponents; p is the gas density; (3 is the temperature expansion 
coefficient; c p is the specific heat of the gas under constant 
pressure; T is the temperature; AT is the temperature differ¬ 
ence; g x and g y are the projections of the free fall acceleration 
vector on the coordinate axes; A,® is the thermal conductivity 
of the porous gas-filled medium without accounting for gas 
movement; for our case A.® = X q . 

The results of the numerical solutions of equation system 
(8) -(11) are represented in the criterial form: 

E k = /(Ra®, h, 8, (p), 


where E k = (A, + A. c )/X q ; cp is the angle of the flat layer (to 
the horizon line); Ra® is the Raleigh filtration number: 


Ra® 


g(3/;Kp 2 c^Ar 

^® 


( 12 ) 


Based on the numerical calculation results, it was estab¬ 
lished that in a horizontal flat heat-insulating layer (cp = 0), 
the effect of convection on the thermal conductivity starts to 
be significant at Ra® > 40. An analysis of the operating tem¬ 
perature conditions in industrial heat-insulating structures 
and the filtration properties of the available fibrous materials 
demonstrated that the maximum value of the Raleigh filtra¬ 
tion number for flat fibrous insulating layers is 30. There¬ 
fore, the effect of convection on horizontally arranged insu¬ 
lation structures can be neglected. For vertically arranged 
structures with h/8 >10 (which condition is virtually always 


where q = r p /r in : r is the pipeline radius (r jn is the insulator 
surface radius). 

In determining Ra® for a circular interlayer using expres¬ 
sion (13), (r in - r ) should be substituted for h. 

Having calculated the values of E k based on equations 
(13) and (14), the convective heat transfer is determined 
from the formula: 


^c = \(E k - 1). (15) 

In order to calculate the convective component of the ef¬ 
fective heat conductivity of fibrous insulation, one should 
know the permeability coefficient K for the particular mate¬ 
rial, which makes up part of the expression determining the 
number Ra® (12). 

The expression for the permeability coefficient calcula¬ 
tion is found used the criterion equation describing the mi¬ 
gration of liquid or gas in a porous medium both in laminar 
and turbulent modes [15]: 

s _ 64 1 

f, v,-1- rj- , (16) 

Re 1.14-21g(l- 0.625m 1 ' 4 ) 


where f m is the resistance coefficient of the porous medium, 
which is part of the calculation formula: 


A P=fm 



HS, 


2 m 


(17) 


where A p is the pressure loss in filtration through a porous 
layer of thickness H\ p is the moving medium density; v f is 
the filtration rate; S is the specific surface area of the po¬ 
rous medium; Re is the Reynolds number for the porous me¬ 
dium; Re = 4mv f /(vS ), where v is the kinematic viscosity. 

In accordance with expression (7), to calculate the pres¬ 
sure loss in the laminar mode, the following formula is used: 

v r |xH 

Ap = ^—. (18) 

Js. 


By equating the right parts of equations (17) and (18) 
and using only the right term in the right part of equation 
(16) for determination of f , we obtain a formula for the cal¬ 
culation of the permeability coefficient for a porous medium: 


R = 



(19) 
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Fig. 2. Thermal conductivity of fiberglass insulating material with 
the mean diameter 5 x 10 " 6 m and the volume weight 60 kg/m 3 ver¬ 
sus temperature: 1, 2, and 3 ) experimental data from [10], [11], and 
[4], respectively; 4) calculation results. 


The specific surface area of fibrous insulation is deter¬ 
mined by the formula: 


_ 4(1-m-n/ p f ) 
sp d c 


( 20 ) 


By substituting expression (20) in formula (19), we find 
the expression for determination of the permeability coeffi¬ 
cient of fibrous insulation: 


K = 


4(1 -m -n / p f ) 


( 21 ) 


where n is the content of nonfibrous inclusions in the fibrous 
insulation, kg/m 3 ; d { is the fiber diameter; in is the fibrous in¬ 
sulation porosity determined by the relationship: 

™ = (pf-p m )/p f , ( 22 ) 

where p f is the fiber density; p m is the fibrous material 
density. 


Using formulae (1), (4) - (6), (12), (15), (21), and (22), it 
is possible to determine the effective thermal conductivity of 
the porous fibrous insulation for any temperature, which is 
corroborated by good convergence between the estimated 
and the experimental data (Fig. 2). 

Thus, a model for estimating heat-shielding properties of 
fibrous insulation made of gas-filled silicate materials is ob¬ 
tained which is suitable for solving practical and engineering 
problems. 
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